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1.0  INTRODUCTION 

The  space  transportation  systems  of  the  future  will  require  new  and  unique 
propulsion  systems.  This  will  especially  be  true  for  the  engines  that  will  be  on  orbit 
transfer  vehicles,  which  will  remain  semi-permanently  in  space.  The  rocket  motors 
on  such  vehicles  will  likely  have  very  high  chamber  pressures,  high  values  of  Isp, 
will  be  throttleable  and  will  have  the  ability  to  start  or  stop  frequently. 

The  orbit  transfer  vehicle  engines  will  also  be  required  to  have  a  nozzle  with  a 
very  large  area  ratio  to  maximize  the  thrust  for  given  chamber  pressure.  The  thrust 
can  be  conveniently  written  in  terms  of  Cpxi  the  thrust  coefficient,  by  the  following 
relation 

T  =  CnPcAt 

where  pc  is  the  chamber  pressure  and  At  is  the  nozzle  throat  area.  The  thrust 
coefficient  is  a  ma.Yimnra  when  the  nozzle  exit  pressure  is  equal  to  the  pressure  of 
the  surrounding  atmosphere  and  when  the  nozzle  is  lossless.  For  lossless  operation, 
Cp  —  cy*,  where  Cf  is  the  thrust  coefficient  for  isentropic  nozzle  flow.  When  the 
external  pressure  is  essentially  zero,  which  would  be  the  case  for  an  orbit  transfer 
vehicle,  then  Cf  reaches  an  asymptotic  limit  when  the  area  ratio,  e  =  Ae/At, 
becomes  very  large  (Ae  is  the  nozzle  exit  area).  Cf  is  seen  to  increase  by  8%  as 
c  is  increased  from  50  to  1000  for  7  =  1.2,  for  example  (Fig.  1).  Thus,  there  are 
advantages  for  increasing  the  area  ratio  if  the  losses  can  be  demonstrated  to  be  low. 

High  expansion  ratio  nozzles  will  create  unique  design  problems  which  have 
not  previously  been  encountered.  The  flow  in  the  nozzle  will  not  be  isentropic, 
especially  in  the  boundary  layers  on  the  wall.  The  boundary  layer  will  be  turbulent 
at  the  nozzle  entrance  but  can  be  seen  to  undergo  a  dramatic  change  in  properties 
before  it  reaches  the  nozzle  exit,  even  if  one  ignores  the  effects  upon  it  of  chemical 
reactions  and  chemical  and  vibrational  relaxation.  During  the  nozzle  expansion 
the  boundary  layer  will  undergo  reverse  transition  and  become  laminar.  By  the 
time  it  has  reached  the  exit  the  flow  will  have  expanded  to  such  an  extent  that 
rarefaction  effects  will  be  noticeable  within  the  boundary  layer.  This  set  of  flow 
circumstances  will  present  a  challenge  to  the  rocket  motor  designer  who  will  have 
to  make  a  very  accurate  estimate  of  the  thrust  coefficient  reduction  to  be  expected 
for  any  particular  nozzle  design. 

Because  the  physics  of  such  boundary  layer  flows  are  so  poorly  understood  it 
will  be  important  to  validate  design  codes  for  large  expansion  ratio  nozzles  by  testing 
prototype  nozzles  before  any  such  nozzles  are  constructed.  Most  conventional  flow 
measurement  tools  will  be  incapable  of  examining  the  boundary  layer  near  the  nozzle 
exit  because  the  local  density  will  be  too  low.  Non- intrusive  sensing  techniques 
(Raman  scattering,  laser  fluorescence,  etc)  will  not  be  usable  inside  the  structure. 
Therefore,  new  techniques  of  observation  will  be  required. 

This  report  describes  the  preliminary  development  of  a  device  which  can  be 
used  to  measure  rarefied  flow  fields,  such  as  the  boundary  layer  near  a  nozzle  exit. 
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The  device  is  a  hot-wire  heat  transfer  probe  which  operates  in  free-molecule  flow; 
it  is  a  modification  of  the  hot-wire  probe  developed  by  Gottesdiener  (Ref.  2).  The 
probe  can  be  used  to  measure  the  flow  properties  within  the  boundary  layer  of  a 
large  expansion  ratio  nozzle  or  other  nozzle  when  the  local  flow  density  is  relatively 
low.  Since  the  wire  diameter  is  10“4  indies,  the  present  probe  can  operate  in  any 
flow  where  the  local  mean  free  path  is  greater  than  5  x  10-4  inches.  It  could 
also  prove  useful  for  measuring  the  flow  within  the  rocket  nozzle  plume  and  in  the 
backflow  region,  upstream  of  the  nozzle  lip.  The  latter  region  is  of  great  interest 
for  spacecraft  contamination  Btudies. 

Although  non-intrusive  techniques  are  normally  used  to  measure  the  plume 
and  backflow  regions,  the  hot-wire  probe  has  some  advantages  over  the  use  of  these 
other  methods.  The  hot-wire  probes  are  relatively  inexpensive  and  arrays  of  fixed 
or  movable  probes  could  be  used  in  these  regions.  Normally  the  alignment  problems 
are  so  severe  for  the  non-intrusive  techniques  that  the  nozzle  is  moved  rather  than 
the  measurement  system.  This  technique  could  not  be  used  for  mpagnring  the  flow 
fields  from  very  large  nozzles.  However,  the  hot-wire  probe  can  be  easily  translated, 
both  inside  the  nozzle  and  outside.  In  addition,  the  probe  can  resolve  the  very  high 
gradients  that  exist  at  the  nozzle  Up  (Refs.  3  and  4)  better  than  the  non-intrusive 
methods  and  can  perform  backflow  measurements  at  lower  densities  than  methods 
such  as  the  electron  beam. 

A  hot-wire  probe  measures  flow  properties  by  measuring  the  change  in  resis¬ 
tance  of  a  heated  filament  which  results  from  the  convective  cooUng  caused  by  the 
flow  over  the  filament.  In  free-molecule  flow  the  probe  can  be  used  to  measure 
the  adiabatic  wall  temperature  and  the  convective  heat  transfer  coefficient  of  the 
filament  wall.  Theoretical  expressions  for  free-molecule  flow  can  then  be  used  to 
relate  the  Stanton  number  and  recovery  factor  (non-dimensional  forms  of  the  con¬ 
vective  heat  transfer  coefficient  and  adiabatic  wall  temperature)  to  the  speed  ratio 
(non-dimensional  velocity).  The  adiabatic  wall  temperature  quickly  approaches  its 
hypersonic  limit  and  cannot  be  used  to  determine  the  speed  ratio.  The  convective 
heat  transfer  coefficient  approaches  its  hypersonic  limit  more  slowly  but  depends 
upon  parameters  in  addition  to  the  speed  ratio.  However,  the  speed  ratio  can  be 
accurately  determined  when  the  hot-wire  probe  measurement  is  coupled  with  an  in¬ 
dependent  measurement  of  the  impact  pressure  or  number  density.  Once  the  speed 
ratio  has  been  determined  then  the  static  temperature,  pressure,  density,  velocity 
and  enthalpy  can  be  calculated.  No  calibration  of  the  probe  is  required  when  it  is 
used  in  free-molecule  flow.  However,  the  probe  surface  area,  the  temperature  coef¬ 
ficient  of  resistivity  of  the  wire  material  and  the  thermal  accommodation  coefficient 
of  the  wire  surface  in  the  flow  must  be  known.  In  addition,  effective  values  of  7  and 
R  for  the  flow  must  be  known  to  calculate  the  flow  properties. 

The  hot  wire  must  be  attached  to  support  posts  which  are  considerably  larger 
than  the  wire.  The  conductive  loss  to  the  supports  is  much  greater  than  the  free- 
molecule  convective  heat  transfer  between  the  heated  wire  and  the  flow.  Very  special 
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design  techniques  were  used  to  attempt  to  eliminate  these  support  conduction  losses. 

The  probe  was  used  to  make  measurements  of  the  flow  field  at  the  exit  of  a  small 
low  density  nozzle.  The  nozzle  configuration  and  nominal  test  conditions  are  shown 
in  Fig.  2.  Probe  measurements  were  made  inside  the  nozzle,  at  several  locations 
downstream  of  the  nozzle  exit,  and  in  the  backflow  region,  upstream  of  the  nozzle 
lip.  Measurements  of  the  heat  transfer  coefficient,  normalized  by  the  centerline 
value,  compared  very  favorably  with  normalized  impact  pressure  measurements  and 
electron  beam  number  density  measurements. 

A  probe  of  this  design  had  been  previously  used  to  measure  a  flow  that  was 
at  ambient  temperature  (roughly  300° K).  By  contrast,  great  temperature  extremes 
existed  in  the  flow  facility  in  which  the  present  probe  was  used.  The  gas  stagnation 
temperature  was  500°K  while  the  free-molecule  recovery  temperature  was  about 
580° K,  which  required  heating  the  probe  wire  to  a  maximum  of  625°K.  On  the 
other  hand,  the  nozzle  walls  were  less  them  500°K  and  the  entire  nozzle  and  flow 
were  surrounded  by  a  77°K  cryogenic  pumping  shroud.  Under  these  conditions  the 
radiation  loss  from  the  probe  surfaces  became  much  greater  than  the  wire  support 
conduction  loss.  Attempts  were  made  to  reduce  this  error  but  they  were  not  entirely 
successful.  Suggestions  for  further  reducing  the  radiation  error  are  given  at  the  end 
of  the  report. 

The  probe  used  for  the  tests  reported  upon  in  this  report  had  L/D  =  2280  and 
was  in  free-molecule  flow  (Kn  >  5)  if  n  <  8  x  1016 /cm3.  For  the  flow  field  into 
which  it  was  inserted  the  Knudsen  number  based  on  the  wire  diameter  was  always 
greater  than  80.  Under  these  conditions  the  probe  was  essentially  non- intrusive  to 
the  flow  field. 

The  background  for  the  development  of  the  free-molecule  heat  transfer  probe 
is  described  in  Section  2.0.  The  theoretical  description  of  the  operation  is  given  in 
Section  3.0  while  the  probe  and  accompanying  electronics  are  described  in  Section 
4.0.  The  measurement  results  are  discussed  in  Section  5.0.  Finally,  suggestions  for 
future  improvements  in  the  probe  are  given  in  Section  6.0. 

2.0  FREE  MOLECULE  HEAT  TRANSFER  PROBE  DEVELOPMENT 

The  probe  that  was  designed  was  a  modification  of  designs  that  have  been 
successfully  used  by  other  investigators  and  ideas  from  various  previous  investigators 
were  incorporated  into  the  probe.  The  evolutionary  development  of  the  probe  will 
be  described  in  this  section  and  mention  will  be  made  of  the  improvements  that 
were  made  over  previous  probes.  The  present  probe  was  used  in  a  more  severe 
environment  (higher  T0,  lower  Too)  than  experienced  by  previous  works,  as  explained 
in  the  Introduction. 

The  cross-stream  cylinder  was  one  of  the  first  shapes  to  be  examined  theo¬ 
retically  in  free-molecule  flow  (Ref.  5)  and  then  the  theory  was  applied  to  obtain 
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flow-field  properties  from  the  heat  transfer  measurements  of  a  heated  cross-stream 
cylinder.  Sherman  (Ref.  6)  measured  Taw  for  a  very  long  unheated  wire  that  was 
passed  through  a  low  density  plane  shock  wave  and  used  this  measurement  to  obtain 
the  shock  wave  profile  and  thickness.  Wong  (Ref.  7)  developed  a  more  traditional 
hot  wire  that  was  used  to  make  measurements  of  hc  and  in  rarefied  flows  at  M 
=  4  and  6.  Coupled  with  an  assumption  of  isentropic  flow  he  was  able  to  use  the 
measurements  to  obtain  the  Mach  number  (speed  ratio)  and  stagnation  tempera¬ 
ture.  Theoretical  expressions  were  used  to  make  corrections  for  conduction  losses 
to  the  supports  and  radiation  losses  were  found  to  be  unimportant  [Taw  ~  345°K, 
Ty,  —  Too  <  100°K].  Laurmann  and  Ipsen  (Ref.  8)  used  this  hot  wire  probe  to  make 
low  density  flow  field  measurements  around  two  simple  two-dimensional  shapes  at 
M  =  4  and  6.  They  use  an  8  inch  long  wire  to  get  a  very  large  length-to-diameter 
ratio  to  minimize  the  conduction  error  and  found  that  large  overheat  ratios  were 
required  to  accurately  measure  hc  and  Tam. 

The  heat  conduction  down  the  support  prongs  was  the  major  problem  with  the 
free-molecule  heat  transfer  probe  described  above  and  prevented  its  adoption  as  a 
standard  low  density  measurement  tool.  In  addition  to  this  problem,  Brown  and 
Thomas  (Ref.  9)  pointed  out  that  more  heat  was  conducted  through  the  supports 
when  the  probe  was  in  a  gas  than  when  it  was  in  a  vacuum  (while  keeping  the  av¬ 
erage  probe  temperature  the  same).  Therefore,  the  conduction  had  to  be  addressed 
before  the  probe  could  be  successfully  used  for  flow  measurements.  While  ignoring 
radiation,  Rajasooria  and  Brundin  (Ref.  10)  developed  a  technique  for  directly 
determining  the  conduction  losses  by  calibrating  the  probe  in  a  free  stream  and  in 
a  vacuum.  Their  procedure  also  eliminated  the  need  for  directly  determining  the 
probe  dimensions  or  average  material  thermal  conductivity.  Lord  (Ref.  11)  greatly 
attended  the  technique  for  determining  the  conduction  end  loss  correction,  again 
ignoring  radiation,  and  extended  the  techniques  to  apply  to  fluctuation  measure¬ 
ments  in  supersonic  continuum  flows  (Ref.  12).  Muthoo  and  Brundin  (Ref.  13) 
used  Lord  s  correction  technique  and  coupled  the  hot-wire  measurements  with  those 
of  the  impact  probe  to  obtain  measurements  of  the  speed  ratio,  pressure,  temper¬ 
ature,  stagnation  temperature,  density,  velocity  and  enthalpy  behind  a  sphere  at 
a  Mach  number  of  5.6.  They  resistively  heated  the  support  posts  to  minimize  the 
conduction  loss. 

A  major  advance  was  made  by  Gottesdiener  (Ref.  2)  who  placed  small  bead 
thermistors  at  the  ends  of  the  heated  wire  and  used  the  thermistors  as  nnnlim«ir  heat 
sources  in  a  temperature  control  system  that  essentially  allowed  the  entire  length  of 
wire  to  be  at  a  constant  temperature.  This  directly  eliminated  the  conduction  loss 
and  allowed  the  radiation  corrections,  which  were  applied  to  each  measurement,  to 
be  measured  under  vacuum  conditions.  He  used  the  probe  in  a  free-jet  flow  where 
M  —  2  and  p  =  2  x  10  5  torr.  Coupled  with  electron  beam  density  measurements, 
he  was  able  to  determine  the  flow  speed  ratio  at  the  probe  location.  His  design 


10 


AEDC-TR-B0-22 


has  been  basically  copied  for  the  probe  that  is  described  in  Section  4.0.  However, 
radiation  presented  a  much  greater  problem  than  was  confronted  by  any  previous 
investigators  and  the  probe  was  radiation  shielded  for  measurements  in  the  lowest 
density  flows,  especially  when  making  measurements  outside  the  nozzle  where  = 
77°K. 


3.0  THEORY 

3.1  Free-Molecule  Heat  Transfer  Tb  a  Cross  Stream  Circular  Cylinder 

The  free-molecule  flow  over  an  infinitely  long  circular  cylinder  at  right  angles 
to  the  flow  field  can  be  determined  by  assuming  that  the  cylinder  wall  temperature 
is  constant  over  the  entire  surface,  that  the  thermal  accommodation  coefficient  is 
constant  over  the  surface  and  that  each,  molecule  carries  j/2  kT  units  of  internal 
energy,  where 

■  _  5~3t 

3~  7-1  ' 

The  thermal  accommodation  coefficient,  a,  is  usually  assumed  to  be  the  same  for  all 
modes  of  energy  transfer  since  experimental  evidence  is  lacking  about  the  accom¬ 
modation  for  internal  energy.  The  flow  is  assumed  to  be  uniform  and  in  Maxwellian 
equilibrium.  Then,  the  Stanton  number  and  recovery  factor  are  given  by  (Ref.  14) 
(see  figures  3  and  4). 


7  (3  +  2.>)/,(4)+(l  +  2.2)/,(^) 
r_7  +  i  (i  +  *2)Mlr)  +  s2/i(<f) 

(1) 

(2) 

These  expressions  can  be  simplified  by  defining  the  following  functions  (Ref.  15) 

Z,  =  (j) 

(3) 

=  tA-"'  ['"  (t)+i‘(y)] 

(4) 

or  the  functions 

/(s)  =  Z!(s2+3)  +  Z2(s2+7/2) 

(5) 

<j(s)  =  3(J?i  +  Z2). 

(6) 
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The  functions  /  and  g  are  plotted  in  Fig.  5.  In  terms  of  these  functions 

r  _  27  a2Z\  +  (a2  +  |)Z2  _  67  f(a)  —  g(s) 

7  +  1  s2(Zt  +  Z2)  7  +  1  s2ff(s) 

of_  7  +  1  1  „  >,  .7  +  1  1  1 K-*) 

Si~a  7  4^s(Z,  +  22)-"“  “T  ’ 


Also  (see  Fig.  6), 


Taw  _  f(s) 

T  g(s)  ■ 


Muthoo  and  Brundin  (Ref.  13)  combined  these  equations  with  the  equation  de¬ 
scribing  the  measurement  of  an  impact  probe 


Pi  =  pK(s) 


to  get  the  result 


For  s  >  ^7/2 


'icy/2irTate  _  ( 7  +  1\ 

3 Piy/R  ~a\y-l) 


V7(*)g(-0 

18^(5)  ' 


Thus,  measurements  of  hc>Taa  and  p{  yield  the  speed  ratio,  if  a  is  known.  The 
function  in  Eq.  (11)  is  plotted  in  Fig.  7. 

In  contrast,  Gottesdiener  (Ref.  2)  combined  the  hot-wire  measurements  with 
an  independent  measurement  of  the  density,  using  the  relation 


P\/Taw 


to  obtain  the  speed  ratio.  The  velocity  was  then  determined  from  a  plot  of 


V2RTZ  -  JV  /W  {  } 

versus  F.  These  relations  are  given  in  Figs.  8  and  9. 

Onoe  the  speed  ratio  is  known,  plus  either  pi  or  p,  the  remaining  flow  properties 
follow  from  Eqs.  (9),(10),  the  definition  of  the  speed  ratio  and  assumptions  of  a 
perfect  gas,  namely, 

P  =  pRT  (15) 
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H  =  -~rRT  .  (16) 

7-1 

Certain  asymptotic  forma  of  Eqs.  (1)  and  (2)  are  required  for  the  data  reduction 
or  are  useful  for  understanding  the  probe  operation.  First,  consider  the  static  gas 
approximation  (s  =  0).  An  expansion  for  the  modified  Bessel  functions,  valid  near 
s  =  0,  is  given  by 


J„<z) 


-GO'S 


(l/4z2) 


'0  Air^  +  i  +  1) 

l/4z 


W-(^)  {j,^^  +  r^  +  2)+”} 


Using  the  perfect  gas  expression  for  c«,  i.e., 

1  a 


Cv  = 


7-1 


R 


and  the  definition  of  Si,  then,  as  s  — *■  0,  we  get  to  lowest  order 

hc  _  wMg  +  1)  =  «£^ r+1)  f1+l  s,/16l 

4y?r  8^2  J 

hc "  a  • 


(17) 


(18) 

(19) 


The  last  form  will  be  used  to  determine  a  from  heat  transfer  measurements  in  a 
static  gas. 

/  ,  v  (3  +  2,j)/.(4)  +  (1  +  2»2)/i(4) 

’’"Vl'+V  (l  +  »J)ii(£)+»2Ji(£) 

Also, 
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The  experiments  were  all  performed 

7) 


at  high  values  of  a.  In  this  case, 


a 


oo  (Ref. 


and 


s<=“(:!Ti)5;(1+4lj+^+”)  •  <*) 

Evaluating  ail  properties  at  the  static  gas  conditions  approaching  the  probe  anrl 
using  the  perfect  gas  expression  for  cpt  i.e. 


leads  to 


R 


or 


Also, 


leading  to 


*■ = “vb  (£4)  {'+i  + 

+  A  +  ^  +  +...) 

/W  =  ^(2’,  +  f‘  +  ^  +  3^+') 


(23) 

(24) 


and 


(25) 

(26) 


A  comparison  of  Eqs.  (23)  and  (26)  indicates  that  Taiu  approaches  its  asymptotic 
limit  much  faster  than  hc  does.  Because  the  wall  temperature  is  nearly  constant  for 
s>  2,  the  heat  transfer  coefficient  must  be  measured  at  high  speed  ratios  to  obtain 
an  accurate  measurement  of  the  local  flow  properties.  However,  a  measurement  of 
he  alone  will  not  result  in  the  speed  ratio.  Therefore,  either  pi  must  be  measured 
and  a  determined  from  Eq.  (11)  or  p  must  be  measured  and  s  determined  from  Eq. 
(13).  External  to  the  nozzle  exit  either  could  be  measured,  p  by  the  electron  beam 
technique,  but  inside  the  nozzle  the  impact  pressure  must  be  measured,  along  with 
the  hot-wire  probe  measurements. 
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3.2  Heat  Transfer  Analysis  of  Probe 

Consider  the  steady  state  heat  transfer  from  a  finite  length  wire  mounted  upon 
supports  that  are  maintained  at  temperatures  Tt  by  end  heaters. 


From  an  energy  balance 

Qj  +  Qt  =  Qc  +  Qr  +  Qk  *  (27) 

It  has  been  assumed  that  the  supports  are  large  compared  to  the  wire  so  that  there 
is  negligible  joulean  heating  of  them.  We  can  assume  one-dimensional  heat  transfer 
in  the  thin  cylindrical  wire  since  L/D  >>  1.  Introducing  the  heat  transfer  coefficient 
leads  to  the  equation 

K  =  Q‘  j  ~Qr  (wrtt»/m»K) .  (28) 

Let  the  wire  resistance  be  given  by 


jftui  —  iJa[l  +  Of(7ui  —  7®)] 


where  Ra  is  the  resistance  at  the  reference  temperature  Ta  and  ar  is  the  temperature 
coefficient  of  resistivity.  Several  reference  temperatures  have  been  tried  and  the  only 
one  that  simplifies  the  analysis  is 


7a  —  TJiu) 


where  Tav>  is  constant  at  each  probe  location. 
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Thus, 

—  iJau([l  +  Ctr(Tw  —  Tm,)]  .  (29) 

Lord  (Ref.  11)  has  examined  the  heat  transfer  to  a  hot  wire  probe  and  shown 
that  it  is  possible  to  obtain  experimental  corrections  for  the  heat  losses  due  to 
conduction  to  the  support  posts.  One  of  the  useful  features  of  his  technique  is 
that  the  effective  wire  length  and  diameter  do  not  have  to  be  measured  but  are 
determined  by  calibration.  However,  he  ignored  radiative  heat  transfer  to  the  probe. 
The  application  to  which  the  present  probe  was  used  does  not  allow  for  the  neglect 
of  radiation.  Therefore,  an  attempt  was  made  to  extend  the  analysis  of  Lord  to 
include  radiative  transfer.  This  analysis  is  included  below,  where  it  is  assumed 
initially  that  Qt  =  0. 

For  a  short  length,  dx,  of  the  wire  the  energy  balance  equation  can  be  written 
as  (Refs.  7  and  11) 

jyi  j  r  -I 

~WT  [*•  ■rffj  +  *Dh(T»  -  T™)  +  ~  T*>) 


fir? 

=  -j^[l  +  C',(Tw-Taw)]  .  (30 

It  has  been  assumed  that  the  wire  is  completely  surrounded  by  a  surface  at  a  uniform 
temperature  Too  and  that  all  surfaces  are  gray. 

A  general  solution  to  Eq.  (30)  cannot  be  written  down  without  making  some 
assumptions,  because  of  the  nonlinearity  of  the  equation.  It  will  be  assumed  that 
the  wire  thermal  conductivity  is  constant  along  the  wire  and  that  the  radiation 
term  can  be  linearized.  Because  the  probe  will  be  operated  with  essentially  uniform 
temperature  over  its  length,  the  first  approximation  is  valid.  The  linearization  of 
the  radiation  term  is  not  justified  for  the  conditions  under  which  the  probe  was 
used  but  must  be  done  to  get  an  analytical  solution  to  the  equation.  However,  in 
many  cases  it  is  justified.  The  method  that  was  used  for  attempting  to  reduce  the 
radiation  term  will  be  discussed  in  Sections  4.3  and  5.0. 


Let  6  —  T„  —  rau„  where  9  will  be  assumed  to  be  small  compared  to  Taxu.  Then, 


T*  —  T4  ~  T4 

«E  t»  —  ^aui 


Using  this  approximation  Eq.  30  becomes 


(31) 


<P8 

dx* 


—  AS  =  —B 


(32) 
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4  4Ac 

DKm 

4  ParRaw  1  BaewT^w 
,  *D2LKW  DKw 

(33) 

4 12Raw 
7rD2LK, 

4eretoT*V)  J"  /  Too  A  1 

*  "  DKw  [  \Taw)  J  ' 

(34) 

This  equation  will  be  solved  using  the  following  boundary  conditions 

1. 

T„  =Tt  at  x  =  ±L/2 

(35) 

2. 

dTw 

— — -  =  0  at  x  =0. 
ax 

The  solution  is 

Tw  =  Tav>  +  ^ 

rr,  rry  COshA.1/2® 

Tt  Taw  a)  cosh A1/2L/2  ■ 

(36) 

This  is  the  same  solution  given  by  Lord  (Ref.  11)  except  that  he  ignored  the 
radiation  term. 


The  measured  wire  resistance  is  given  by 

=  J%~ f  P  +  «r(Tw  -  Taw)]dx 

„  „  „  o,rR„„B\  tanl.4‘^1/2) 

+  A  Ita  w  A  J  (A1/iLf2) 

For  an  infinitely  long  wire  the  resistance  is  given  by 


Rw„  =  Raw  + 


OlfRawR 


Note  that 


,  ,  t,  n  4 hc  1  SctJTln  4arcrewT*„ 

A  +  arti  =  G  =  -W7T-  + 


DKU 


DKU 


DK„ 


1  - 


Then 

We  can  write 


RfVca  —  Q  R^Voa  d"  Raw  • 


Of -B  /  4 ctrRaw  \  2  4ar 

~C~  ~  \C*D2LKW )  C 

=  xi2  -  r; 


&^w  Til 


(ft)')  • 


CDKU 


-m 


(37) 


(38) 


(39) 

(40) 


(41) 
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then 

a*.  =  (jjy) +  (tTt)  r 

Thus,  RWaa  versus  is  a  straight  line  with  slope 


and  intercept  (i+y )  Rat 


m0 o 


X 

i  +  r 


(42) 


(43) 


where 


m°°  =  7'DL(hc+4acwT3J 

^  i  _  Q'rP'wem^u,[l  ~  (Tpo/Jam)1] 

1+Y  hc  +  ±otJTlw 


(44) 

(45) 


These  are  the  expressions  given  by  Lord  (Ref.  11)  with  the  addition  of  the  radiation 
terms. 


Thus,  it  is  observed  that  the  technique  of  Lord  for  correcting  for  the  conduction 
to  the  supports  cannot  be  used  unless  the  radiation  term  can  be  neglected  compared 
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to  the  convection.  Consider  a  typical  condition  near  the  exit  of  the  Mach  number 
—  4  low  density  nozzle  used  during  the  measurements: 

gas  -C02,n  =  10 21/m\Taw  =  583“K  .T*,  =  77°K  , V  =  815 m/s. 


Then, 


fee 

4cr  euT*w 


2.26. 


Flrom  Eq.  (23),  fec  diminishes  as  nV  diminishes  and  the  radiation  correction  becomes 
even  greater  as  the  probe  is  moved  downstream  into  the  nozzle  plume. 


This  large  magnitude  of  the  radiation  contribution  invalidates  the  correction 
procedure  of  Lord  (Ref.  11).  In  the  present  case,  it  was  attempted  to  eliminate  the 
conduction  correction  by  placing  thermistors,  used  as  heaters,  at  the  ends  of  the 
wire  and  applying  power  such  that 


Qt  =  Qk  • 


(46) 


Assuming  that  this  is  possible  Eq.  (28)  becomes 


K  = 


Qi-Qr 

Am(Tw  —  Tsui) 


watts/m2oK  . 


(47) 


Radiation  was  a  considerably  greater  problem  for  the  probe  when  used  for  the 
present  tests  than  for  the  tests  of  previous  investigators  because  of  the  high  stagna¬ 
tion  temperature  and  correspondingly  high  adiabatic  wall  temperature  (500°  K  and 
583°K,  respectively)  and  because  the  walls  of  the  chamber  were  cooled  with  liquid 
nitrogen.  A  simple  radiation  heat  transfer  analysis  of  the  probe  would  assume  that 
the  probe  was  surrounded  by  walls  at  a  temperature  Too  and  that  all  walls  were 
gray.  The  rate  of  radiative  heat  transfer  would  then  be  given  by 

Qr  =  ewV(i:-^)  (48) 

where  ew  is  the  emissivity  of  the  probe  (less  than  0.1  for  the  present  probe). 

One  way  to  compare  the  influence  of  radiation  on  the  probe  is  to  compute  the 
ratio  of  the  radiation  to  convection  heat  transfer  rates,  QT/Qc>  given  for  large  speed 
ratios  by  the  ratio  of  Eq.  (48)  to  (23),  i.e. 

Qr  _  _ O’CwCT*  —  T^g) _ 

Q-  ~  (aAj 2n  (^±i)  r»VJb(l  +  1/4 s*){Tw  -  Taw)  ' 
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This  ratio  becomes  large  under  the  following  circumstances:  nV  becomes  small,  Tw 
much  greater  than  Too,  or  Tw  approaches  Taw. 

First,  consider  the  situation  when  the  product  nV  becomes  small.  This  case  is 
shown  in  Fig.  10,  for  particular  values  of  the  other  parameters.  Radiation  is  less 
than  convection  to  the  right  of  the  curve  and  becomes  important  at  a  higher  value 
of  nV  as  Tw,  decreases.  At  the  centerline  of  the  nozzle  for  the  present  tests  (see 
Fig.  2)  nV  =  4.8  x  1024/m2s.  Thus,  at  that  point  Qc  >  Qr.  However,  at  the  nozzle 
lip  the  velocity  was  less  and  the  number  density  less  by  a  factor  of  at  5,  than  at 
the  centerline,  so  that  Qr  >  Qe.  This  also  occurred  in  the  backflow  region,  where 
the  number  density  was  considerably  lower.  The  problem  was  more  severe  in  the 
plume  and  backflow  region  where  Too  —  77°  K  for  a  considerable  Bolid  angle,  whereas 
Too  —  47CTK  within  the  nozzle.  This  problem  is  accentuated  as  Tw  approaches  Taw 
(low  overheat).  Note  that  there  must  be  an  error  in  the  paper  by  Gottesdiener  (Ref. 
2)  who  indicated  that  Qr  =  QC  at  n  —  4  x  1017/m3;  the  present  analysis  indicates 
that  this  occurs  at  much  higher  values  of  number  density  is  also  dependendent  upon 
the  velocity. 

For  a  given  value  of  nV,Qr/Qc  increases  as  Tw  approaches  Taw,  i.e.,  as  the 
overheat  ratio  decreases.  This  effect,  also  noticed  by  Laurmann  and  Ipsen  (Ref.  8), 
is  illustrated  in  Fig.  11,  which  gives  Qr/Qc  as  a  function  of  Tw  for  given  values 
of  the  other  parameters.  It  will  be  shown  in  Section  4.3  that  the  probe  must  be 
operated  at  a  number  of  overheat  ratios  and  the  Rw  versus  (Qj  —  Qr)  curve  must  be 
extrapolated  to  zero  overheat.  The  radiation  correction  and,  correspondingly  the 
uncertainty  in  the  measurement,  will  increase  as  the  overheat  ratio  dprrrews 

The  above  analysis  has  considered  the  radiation  only  from  the  heated  filament 
However,  potentially  the  major  source  of  radiation  loss  is  from  the  surface  of  the 
thermistors.  The  thermistors  and  thermocouples  were  attached  to  the  support  posts 
with  a  high  temperature  casting  resin  which  created  a  black  spherical  blob  at  the 
end  of  the  wire  which  had  a  diameter  that  was  estimated  to  be  less  than  0.020  inches 
in  diameter  (see  Fig.  14).  The  surface  area  of  the  resin  spheres  was  up  to  35  timpg 
larger  than  that  of  the  wire  and  could  correspondingly  radiate  a  couple  of  orders 
of  magnitude  more  power  than  the  wire.  Evidence  that  this  could  actually  be  the 
case  is  given  in  Section  5.0  and  suggestions  for  significantly  reducing  this  loss  will 
be  given  in  Section  6.0. 

4.0  PROBE  DESCRIPTION  AND  OPERATION 
4.1  Probe  Construction  and  Calibration 

The  probe  construction  is  shown  schematically  in  Figs.  12  and  13  and  a  pro¬ 
totype  of  the  probe  is  photographed  in  Fig.  14.  Two  probes,  having  nominal  wire 
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lengths  of  1/8  and  1/4  inches,  were  constructed.  The  heated  element  consisted  of 
a  100  x  10-6  inch  diameter  Pt-10%  Rh  wire.  This  wire  was  manufactured  by  the 
Wollaston  process  and  had  a  silver  sheath  that  was  removed  by  a  50%  nitric  add 
solution.  The  wire  was  attached  to  0.003  inch  diameter  constantan  posts  which 
were  in  turn  mounted  to  0.020  inch  diameter  constantan  supports.  Normally  the 
wires  would  be  attached  to  nickel  supports  but  constantan  has  one-fourth  the  ther¬ 
mal  conductivity  of  nickel  and  its  use  minimized  the  conduction  loss  from  the  ends 
of  the  wire.  The  constantan  posts  were  silver  soldered  to  the  .020  inch  diameter 
constantan  supports  with  a  eutectic  silver  solder  and  the  Pt-10%  Rh  wire  was  spot 
welded  to  the  0.003  inch  diameter  constantan  posts  using  the  spot  welder  shown  in 
Fig.  15. 

Following  Gottesdiener  (Ref.  2),  thermistors  were  placed  at  the  ends  of  the 
wire  in  an  attempt  to  eliminate  the  conduction  losses  by  heating  the  ends  of  the 
wire  an  amount  so  that  Qt  =  Qk-  Thermistors  are  useful  heating  elements  to  use  in 
proportional  temperature  controllers  because  of  their  rapid  variation  of  resistance 
with  temperature.  The  thermistors  used  were  VECO  0.013  inch  diameter  bead 
thermistors  with  a  resistance  of  1000  ohms  and  dissipation  constant  of  O.10mu>/°C. 
They  could  deliver  a  maximum  power  of  0.012  watts.  The  temperature  variation  of 
the  thermistor  resistance  is  shown  in  Fig.  16.  The  end  temperatures  were  sensed 
with  0.002  inch  diameter  chromel-alumel  thermocouples.  The  thermocouple  wires 
were  insulated  and  enclosed  in  a  0.010  inch  outside  diameter  stainless  steel  sheath. 
The  thermocouples  and  thermistors  were  attached  to  the  constantan  posts  with 
Stycast  2762  FT,  which  is  a  high  temperature  casting  resin  that  possesses  a  high 
thermal  conductivity.  Before  attaching  the  thermocouples  and  thermistors,  the 
constantan  posts  were  coated  with  Cal  Custom  #8028  white  high  temperature  paint. 
The  Pt-10%  Rh  wires  were  attached  after  the  thermocouples  and  thermistors  were 
in  place. 

The  constantan  supports  were  glued  into  a  ceramic  support  that  consisted  of 
a  cluster  of  multi-holed  ceramic  tubes.  The  ceramic  tubes,  in  turn,  were  mounted 
inside  a  stainless  steel  tube.  The  thermocouple  and  thermistor  lead  wires  were 
led  through  the  remaining  holes  in  the  ceramic  tubes.  See  Figs.  12  and  13  for 
construction  details. 

The  stainless  steel  tube  was  supported  on  a  traverse  inside  tbe  vacuum  chamber 
that  allowed  it  to  be  rotated  and  translated  with  respect  to  the  plane  of  symmetry 
of  the  nozzle.  The  short  length  of  the  wires  allowed  the  probes  to  be  positioned 
to  within  several  thousandths  of  an  inch  of  the  nozzle  wall.  This  was  more  than 
adequate  to  resolve  the  properties  through  the  nozzle  wall  boundary  layer  and  their 
diameter  was  small  enough  to  resolve  the  high  flow  gradients  in  the  nozzle  lip  region. 

The  probes  were  constructed  and  calibrated  by  personnel  of  Calspan  Field 
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Services,  Inc.,  AEDC  Division,  at  Arnold  Air  Force  Station.  It  is  very  crucial 
that  the  probe  resistance  and  temperature  coefficient  of  resistivity  be  known  very 
accurately  (to  0.1%)  to  accurately  use  the  probe  to  measure  flow  properties.  The 
equipment  used  to  calibrate  the  probe  is  shown  in  Fig.  17.  The  probe  and  an  RTD 
probe  were  placed  in  a  large  steel  block  which  was  placed  in  an  oven.  The  oven 
temperature  was  increased  to  greater  than  500°F  and  the  oven  power  was  turned 
off.  The  RTD  probe  was  connected  to  a  Fluke  2180A  RTD  digital  thermometer, 
whose  output  was  fed  to  the  HP  computer.  Data  were  taken  when  the  oven  cooled 
to  500°  and  at  approximately  each  10°F  decreased  until  the  temperature  reached 
80°F. 

The  probe  wire  and  each  thermistor  (TM1  and  TM2)  were  attached  to  the 
voltage  divider  circuit  shown  in  Fig.  17,  where  the  value  of  the  divider  resistor  is 
indicated.  The  two  thermocouples  (TCI  and  TC2)  were  connected  to  a  Kaye  ice 
point  reference.  All  voltages  were  measured  with  a  Fluke  8502A  voltmeter  after 
passing  through  a  Fluke  2204A  scanner.  The  four  wire  technique  was  used  to 
measure  all  the  voltages.  Input  voltages  of  approximately  2, 4,6, 8,  and  10  volts  were 
applied  to  the  voltage  divider  circuits.  Thirty-two  readings  were  internally  averaged 
at  each  reading.  The  procedure  at  each  oven  temperature  was  as  follows: 

1.  Read  the  oven  temperature. 

2.  Read  the  voltage  from  TCI  and  TC2. 

3.  Read  the  oven  temperature. 

4.  Read  the  wire  potential  for  the  five  values  of  Ei. 

5.  Read  the  oven  temperature. 

6.  Read  the  TMl  potential  for  the  five  values  of  £;. 

7.  Read  the  oven  temperature. 

8.  Read  the  TM2  potential  for  the  five  values  of  E,. 

9.  Read  the  oven  temperature. 

The  calibration  results  are  shown  in  Figs.  13  through  21.  The  probe  wire 
resistances  are  shown  in  Figs.  18  and  20  for  the  1/8-inch  and  1/4-inch  probes, 
respectively,  and  the  thermal  coefficient  of  resistivity  was  determined  from  a  least 
square  fit  through  the  measured  points.  The  two  determined  values  of  ar  agree 
closely.  The  measured  thermistor  resistances  are  shown  in  Figs.  19  and  21.  The 
measured  thermocouple  outputs  corresponded  closely  to  the  standard  tabulated 
values. 

At  very  low  values  of  the  product  nV  the  radiation  heat  transfer  from  the 
heated  wire  is  much  larger  than  the  convection.  When  the  probe  was  used  in  these 
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circumstances  it  had  been  planned  to  surround  it  by  a  double  cylindrical  radiation 
shield  shown  in  Fig.  22.  The  inner  Bhield  was  to  consist  of  two  layers  of  copper  wire 
(30  mesh,  37%  open  area)  which  would  have  been  laid  at  45“  to  one  another.  This 
shield  would  have  contained  an  insulated  ni chrome  wire  which  would  have  heated 
it  to  the  wire  temperature.  If  the  shape  factor  of  this  shield,  from  the  probe,  had 
been  1.0,  then  heating  the  shield  to  Tw  would  have  eliminated  the  radiation  error. 
However,  the  probe  would  no  longer  have  been  in  f roc-molecule  flow.  Therefore, 
the  wire  mesh  was  to  have  been  used;  the  molecules  which  struck  the  mesh  would 
have  been  scattered  in  many  directions  with  little  probability  of  returning  to  the 
wire,  and  those  that  proceeded  through  the  mesh  would  been  primped  by  the  outer 
cylinder,  which  was  to  be  cooled  with  liquid  nitrogen  (CO2  is  effectively  pumped 
at  77“K).  This  shielded  arrangement  would  have  allowed  the  probe  to  have  been 
used  in  regimes  never  reached  by  previous  investigators  but  it  was  never  constructed. 
Therefore,  it  is  not  known  how  effective  it  might  have  been  in  reducing  the  radiation 
loss.  Instead,  a  simple  shield  consisting  of  a  single  helical  coil  of  0.3  inch  diameter 
Nichrome  wire  was  placed  around  the  probe.  A  thermocouple  that  had  been  silver 
soldered  to  the  Nichrome  wire  was  used  to  measure  the  shield  temperature  and 
control  the  current  through  the  wire.  This  shield  was  not  very  effective  in  reducing 
the  radiation  loss.  Experience  with  this  simple  radiation  shield  will  be  discussed  in 
Section  5.0. 

Seb&n  (Ref.  16)  measured  the  emisaivity  of  Pt-13%  Rh  for  temperatures  greater 
than  700°K  (Fig.  23).  He  found  that  the  emissivity  greatly  increased  if  the  Rhodium 
oxidized  (x  data  points).  Care  was  taken  to  avoid  wire  oxidation  during  the  manu¬ 
facture  and  testing  of  the  probe  and  the  wire  always  seemed  to  maintain  its  shiny 
appearance.  At  the  operating  temperatures  used  in  the  experiments,  the  data  in 
Fig.  23  indicates  that  eu  <  0.1.  At  any  temperature  the  emissivity  can  be  approxi¬ 
mated  by  the  equation 

ew  =  0.57Qy/pwTw  -  0.179(pt.T.)  +  0 M(pwTw)3'2  (51) 

where  pw  is  the  wire  resistivity,  in  units  of  ohm  cm,  and  Tw  has  the  units  of  “K. 
The  resistivity  is  given  by 


Pw  =  Po[l  +  <*r(Tu>  —  273)].  (52) 

For  Pt-10%  Rh,  po  =  18.4  x  10-6  ohm,  cm  and  ar  =  0.0017/“K.  At  Tw  =  625°K, 
this  formula  gives  ew  =  0.076.  A  conservative  value  of  0.1  has  been  used  in  all  the 
heat  transfer  calculations. 

The  heat  transfer  from  the  wire  also  depends  on  the  thermal  accommodation 
coefficient  of  the  wire.  The  accommodation  of  both  the  translational  and  internal 
energies  must  be  considered.  Very  little  evidence  exists  about  the  accommodation 
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of  internal  energy  but  that  which  does  exist  indicates  that  the  accommodation  co¬ 
efficients  for  CO2  for  translational  and  rotational  energy  are  both  equal  to  1.0  while 
the  accommodation  coefficient  for  vibrational  energy  is  small,  so  that  the  overall 
accommodation  coefficient  is  0.81  (Ref.  17).  The  possibility  existed  for  nonequilib¬ 
rium  to  occur  in  the  exhaust  of  the  nozzle,  especially  in  the  boundary  layer  and  the 
lip  and  backflow  regions.  Since  the  lowest  vibrational  temperature  of  CO2  is  961°K 
and  the  stagnation  temperature  is  only  500°K,  there  will  be  very  little  excitation  of 
any  of  the  vibrational  modes,  especially  at  the  nozzle  exit  where  the  temperature 
is  about  120s K,  and  only  rotational  nonequilibrium  needs  to  be  considered.  As¬ 
suming  equal  values  of  the  translational  and  rotational  accommodation  coefficients, 
Sherman  (Ref.  6)  derived  the  free-molecule  result 


_  /(3)  -  1  \ 
Tt  ",W  +  3U  7 


(53) 


which  gives  the  same  result  as  Eq.  (9)  when  the  translational  static  temperature 
Tt,  equals  the  rotational  temperature,  Tr.  Using  this  result,  he  calculated  the 
difference  in  Tauj  downstream  of  a  shock  wave  in  a  diatomic  gas,  a-sauming  that 
Tr,  remained  at  the  upstream  temperature,  and  compared  it  to  the  result  that 
would  be  obtained  for  equilibrium  flow.  Even  at  an  upstream  Mach  number  of  4, 
the  difference  in  wall  temperature  was  only  about  2%.  Thus,  even  in  such  a  severe 
nonequilibrium  flow  as  a  normal  shock  wave,  the  effect  of  rotational  non-equilibrium 
on  the  wire  temperature  can  be  neglected  as  long  as  the  accommodation  coefficients 
for  translational  and  rotational  energy  are  equal. 

A  greater  error  than  that  caused  by  the  nonequilibrium  flow  is  caused  by  flow 
gradients  along  the  wire  length.  Traditionally  free-molecule  wire  probes  have  been 
long  (up  to  8  inches)  to  minimize  the  conduction  loss.  The  present  probes  were 
considerably  shorter  and  therefore  could  be  used  in  regions  of  larger  temperature 
gradients. 


The  following  heat  transfer  rate  formulas  were  used  to  assist  in  the  design  of  the 
electronics  and  the  experiments.  They  apply  to  the  1/4-inch  wire  for  the  conditions: 
gas  -  CO2 ,  V  -  815 m/s,  n  =  5.85  x  1021]/m3,s  =  3.25,  a  =  0.9,  ew  =  0.1. 


Qc  =  2.94  x  lO-6^  -  Taw]  watts  (53) 

Qr  =  2.88  x  10-lfl  [7£  -  7^]  watts.  (54) 

An  upper  estimate  of  the  radiation  from  the  resin  surrounding  the  two  thermistors 
is  given  by  the  following  formula. 


Qr  (thermistors)  «  9.1  x  10-H[T4  -  7£>]  watts  (55) 
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where  T  is  the  unknown  surface  temperature  of  the  resin.  This  formula  assumes 
that  the  diameter  of  the  black  resin  is  0.020  inches  and  the  emissivity  of  the  surface 
is  1.0. 


4.2  ELECTRONICS 

4.2.1  Description  of  Electronics  and  Operational  Procedures. 


Two  separate  electronic  circuits  were  required  to  operate  the  probe,  the  first 
to  maintain  the  wire  at  a  prescribed  resistance  and  the  second  to  maintain  the  ends 
of  the  wire  at  the  temperature  which  corresponds  to  the  prescribed  wire  resistance. 
The  circuit  schematics  are  given  in  Figs.  24  to  26  and  photographs  of  the  completed 
instruments  are  given  in  Figs.  27  through  29. 

The  wire  was  maintained  at  a  constant  mean  temperature  with  an  open  loop 
constant- current  anemometer  circuit  and  was  used  much  as  a  constant-current 
anemometer  is  used  to  measure  flow  temperature.  The  circuit  was  designed  to 
meet  the  special  needs  of  the  present  probes  which  included  greatly  increased  probe 
resistance  (100  to  300  ohms,  compared  to  usual  hot  wire  probe  resistances  of  5  to 
10  ohms)  and  the  need  to  operate  at  very  small  overheat  ratios  which  require  small 
amounts  of  power  (to  10-6  watts).  A  20:1  bridge  was  used  to  increase  the  sensi¬ 
tivity  of  the  measurement  and  the  probe  was  placed  in  one  leg  of  that  bridge  (Fig. 
24).  Initially  a  closed  loop  system  was  designed  but  the  bridge  imbalance  required 
to  obtain  stable  feedback  behavior  prevented  that  technique  from  yielding  accurate 
measurements  of  the  probe  resistance. 

The  hot  wire  control  is  operated  in  the  following  fashion. 

1.  Ground  both  inputs  to  the  current  source  amplifier  and  balance  the  am¬ 
plifier. 

2.  Place  a  calibration  resistor  in  place  of  the  hot  wire  [100.04  ±.0U1  resistor 
for  1/8  inch  probe;  200.55  ±0.1J1  resistor  for  1/4  inch  probe]  and  set 
the  thumbwheel  resistor  to  the  value  to  balance  the  bridge  [200  and  401, 
respectively].  Offset  the  current  amplifier  by  0.100  volts  and  adjust  the 
bridge  imbalance  measuring  amplifier  until  its  output  is  zero,  as  read  by 
M4. 

3.  Replace  the  calibration  resistor  with  the  hot  wire  probe. 

4.  A  desired  probe  resistance  can  be  obtained  by  dialing  the  appropriate  value 
on  the  thumbwheel  resistor  and  changing  the  voltage  offset  to  the  current 
source  until  the  bridge  is  balanced  (meter  reading  M4  is  zero). 
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The  current  limit  is  set  by  placing  the  meter  at  position  M4,  increasing  the  offset 
voltage  on  the  current  source  until  the  maximum  desired  voltage  is  just  exceeded, 
and  then  adjusting  the  potentiometer  on  the  FET  until  the  voltage  just  begins  to 
decrease.  Generally  the  bridge  top  voltage  was  adjusted  to  limit  the  probe  current 
to  about  0.6  ma,  although  the  probes  were  able  to  handle  20  ma  in  air  without 
being  destroyed. 

Because  of  the  large  conduction  heat  loss  from  the  wire  to  the  support  posts, 
a  large  temperature  gradient  will  exist  in  the  wire  until  power  is  applied  to  the 
thermistors.  The  purpose  of  the  proportional  temperature  controllers  is  to  apply 
sufficient  power  to  the  thermistors  to  bring  the  wire  ends  to  the  temperature  which 
corresponds  to  the  mean  wire  resistance,  as  set  by  the  wire  control  (see  eq.  (58)). 
This  is  accomplished  by  the  following  steps  (see  Fig.  26). 

1.  Balance  the  input  to  the  two  amplifiers.  In  switch  positions  Ml  and  M2  the 
input  to  the  amplifiers  are  grounded  and  the  outputs  are  open.  Position  Ml 
is  used  to  balance  the  first  stage  and  M2  the  second  Btage  of  amplification. 

2.  Measure  the  output  of  the  thermocouple  transducer  in  switch  position  M4 
with  the  loop  still  open  find  the  amplifier  inputs  grounded.  Set  the  offset 
to  the  first  stage  (10  turn  potentiometer)  to  a  voltage  just  slightly  greater 
than  the  thermocouple  transducer  voltage. 

3.  Close  the  loop  and  set  the  desired  temperature  by  adjusting  the  offset  to 
give  the  desired  temperature  (see  Eq.  (62)). 

The  maximum  power  that  the  thermistors  can  dissipate  is  approximately  0.012 
watts.  The  power  is  limited  by  attaching  a  resistor  of  appropriate  size  at  the  ther¬ 
mistor  input  (the  resistor  magnitude  depends  on  the  expected  operating  tempera¬ 
ture),  increasing  the  offset  voltage  on  the  first  amplifier  stage  until  the  maximum 
possible  voltage  is  achieved  at  the  upper  end  of  the  resistor  and  then  adjusting  the 
potentiometer  on  the  FET  transistor  until  the  voltage  begins  to  drop. 


4.2.2  Transfer  Function  and  Accuracy  of  The  Electronics. 


If  the  resistors  in  the  bridge  were  perfectly  balanced  then  the  bridge  output  as 
measured  by  the  voltage  across  the  bridge  (M3)  would  be  zero  when  the  thumbwheel 
resistor  was  precisely  20  times  the  wire  resistance.  This  could  be  determined  by 
placing  a  small  current  through  the  bridge  from  the  current  source  and  measuring 
the  voltage  M3.  However,  no  attempt  was  made  to  perfectly  match  the  resistors 
and  instead  the  transfer  function  of  the  bridge  was  measured.  First,  the  bridge  was 
balanced  as  outlined  in  steps  1  and  2  in  the  previous  section.  Then  the  calibra¬ 
tion  resistor  was  replaced  with  a  variable  resistor  called  Dekabox,  manufactured  by 
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Electro-Measurements,  Inc.  The  resistance  values  from  the  Dekabox  were  measured 
with  a  Fluke  8040A  multimeter  and  found  to  have  an  accuracy  of  0.3  to  -0. 117  ,  as 
the  total  resistance  varied  from  100  to  38011,  respectively.  The  offset  voltage  on 
the  current  source  ( E 2)  was  placed  at  0.800,  1.200  and  1.600  volts  and  although 
small  differences  were  seen  as  the  offset  voltage  varied  they  were  much  smaller  than 
the  overall  imbalance  condition  of  the  bridge  and  all  measurements  were  averaged 
together. 

The  transfer  function  was  determined  by  setting  the  Dekabox  resistance  value 
(range  100  to  20011  for  the  1/8  inch  probe  and  200  to  38011  for  the  1/4  inch  probe) 
and  determining  the  value  of  the  thumbwheel  resistor  that  balanced  the  bridge.  The 
measurements  are  shown  in  Figs.  30  and  31  and  indicate  that  a  linear  regression 
line  through  the  measurements  was  adequate.  The  results  are  as  follows: 

A)  100.0411  calibration  resistor 

Ri  (input)  =  1.38  4-  0.96786J?m  (measured  on  thumbwheel)  (56) 

Mean  of  deviation:  0.1611 
Standard  deviation:  0.1011 

B)  200.5511  calibration  resistor. 

fZ,  (input)  =  4.21  +  0.95524 Rm  (measured  on  thumbwheel)  (57) 

Mean  of  deviation:  0.2611 
Standard  deviation:  0.31H 

The  resolution  of  the  variable  thumbwheel  resistor  is  0.5H  (10/20).  Therefore,  both 
regression  lines  are  within  the  resolution  of  the  instrument  and  the  basic  resistance 
measuring  accuracy  of  the  instrument  is  comparable  to  the  instrument  resolution. 
For  those  run  conditions  when  E?  was  less  than  0.800,  the  following  equation  was 
used  to  determine  the  input  resistance. 

Ri  (input)  —  C  X  Rf „  (measured  on  thumbwheel)  (58) 


Ej  (volts) 

c 

0.200 

0.9890 

0.300 

0.9845 

0.400 

0.9810 

0.500 

0.9790 

0.600 

0.9776 

0.700 

0.9768 
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The  probe  resistance  is  given  by  the  equation 


Rw  =  i?2o[l  +  <*w(Tw  —  293)]  (59) 

where  the  reference  state  is  taken  to  be  20°  C.  The  uncertainty  in  temperature  is 
related  to  the  uncertainty  in  resistance  by  the  relation 


AT  = 


A  R 
Q  20^20 


(60) 


Using  the  conservative  estimate  A R  —  Q.5J1  and  the  values  of  i2jo  and  <*20  obtained 
from  the  probe  calibration  (Figs.  18  and  20)  then 


AT -I  26°K  C1/8  Probe) 

\  1.5®K  (1/4  inch  probe)  . 

Now  consider  the  accuracy  of  the  proportional  temperature  controllers.  The  basic 
Analogue  Devices  chip  AD595CD  has  an  accuracy  of  ±1°K  and  this  will  be  the 
accuracy  of  the  controllers  if  other  uncertainties  are  not  added  to  the  system.  First, 
the  thermocouple  outputs  were  measured  for  both  thermocouples  on  each  probe 
during  the  probe  calibration  and  their  outputs  were  found  to  agree  closely  with  the 
standard  tabulated  values.  Therefore,  the  tabulated  thermocouple  voltages  were 
assumed  to  accurately  represent  the  signal  from  the  thermocouples  on  the  probes 
and  were  used  to  determine  the  temperature,  without  correction. 


The  transfer  function  of  the  device  AD595CD,  as  given  by  the  manufacturer 
is, 

E  =  247.3 [(mu)  + 11  x  10-8]  volts  (61) 

where  E  is  the  output  voltage  of  the  device  for  a  Type  K  input  voltage  (mu).  E  was 
measured  to  ±lmu.  This  equation  was  solved  for  the  (mu)  input  from  the  thermo¬ 
couple  and  the  standard  tables  were  used  to  determine  the  temperature.  The  device 
accuracy  of  ±1°C  represented  the  overall  accuracy  with  which  the  temperature  of 
the  ends  of  the  wire  could  be  set,  which  at  times  was  not  accurate  enough  because 
of  the  low  overheat  values  that  were  used,  especially  when  there  was  a  very  small 
amount  of  convective  cooling  (see  Section  5.0). 


4.3  Probe  Operation 


The  free-molecule  heat  transfer  probe  measures  the  convective  heat  transfer 
coefficient  and  the  adiabatic  wall  temperature  of  the  wire.  As  previously  discussed, 
kinetic  theory  can  yield  theoretical  expressions  relating  these  quantities  (i.e.,  he  and 
Tav,)  to  the  properties  of  the  flowing  fluid.  A  measurement  of  Taui  could  be  used  to 
determine  the  speed  ratio.  However,  as  seen  from  Eq.  (26)  Taw  is  very  insensitive 
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to  the  speed  ratio  for  large  speed  ratios  and  cannot  be  usefully  used  to  determine 
the  speed  ratio.  Likewise,  the  speed  ratio  cannot  be  directly  determined  from  a 
measurement  of  hc  because  the  convective  heat  transfer  coefficient  depends  upon 
flow  properties  in  addition  to  the  speed  ratio,  as  seen  in  Eqs.  (23)  and  (24).  There¬ 
fore,  to  obtain  an  expression  from  which  the  speed  ratio  can  be  directly  determined 
with  accuracy  at  a  high  speed  ratio,  an  independent  measurement  must  be  made 
of  another  quantity,  either  the  impact  pressure  [Eq.  (11)]  or  the  density  [Eq.  (13)]. 
The  procedures  for  obtaining  the  remaining  flow  properties  once  the  speed  ratio  has 
been  determined  were  outlined  in  Section  3.1. 

The  convective  heat  transfer  coefficient  and  the  adiabatic  wall  temperature  of 
the  probe  are  determined  at  each  probe  location  as  follows. 

A  desired  wire  temperature  is  obtained  by  setting  the  bridge  to  the  resistance 
corresponding  to  the  wire  temperature,  according  to  Eq.  (59),  namely 


Rw  =  T20(l  +  ar(Tw  -  293)]. 


(59) 


Tw  should  be  greater  than  Taw  but  the  ratio  Tw/Taw  had  to  be  less  than  about  1.07 
because  the  thermistors  and  glues  could  only  operate  to  about  625°K.  Next,  power 
is  applied  to  the  thermistor  heaters  to  bring  the  ends  of  the  wire  to  the  tempera¬ 
ture  Tw.  This  step  was  not  possible  for  the  experiments  that  were  performed,  as 
explained  in  Section  5.0.  At  this  condition  Qt  =  Qt  and  Eq.  (28)  becomes 


h  -  Qi 
c  AS(TW  -  Taw) 

Rom  the  hot  wire  circuit  diagram  (Fig.  24), 


(62) 


(63) 


where  72^  is  determined  from  the  regression  Eqs.  (56)  or  (57).  Qr  can  be  determined 
by  repeating  the  same  measurement  of  Qj  under  static  conditions  at  a  very  low 
vacuum  (10-5  torr  or  less).  When  the  two  measurements  of  Qj  are  of  comparable 
magnitude  then  the  probe  must  be  surrounded  by  the  radiation  shield  and  the 
experiments  repeated.  By  replacing  (Tw  -  Tato)  by  Eq.  (59)  we  find  that 


and 


(72w  72au ,)  — 


#20^20 

hcA, 


(Qi-Qr) 


(64) 


dRu,  R20O20 

d{Qj~Qr)~~^AT  ■ 


(65) 
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Therefore,  hc  is  obtained  from  the  slope  and  RaW  is  given  by  the  intercept  of  the 
curve  of  Rw  versus  (Qj  —  Qr).  Taw  is  then  determined  from  Eq.  (59).  At  least  five 
measurements  must  be  made,  at  different  values  of  Rw,  both  with  and  without  flow 
(i.e.,  vacuum),  to  obtain  an  accurate  measurement  of  the  slope  of  the  curve.  This 
procedure  must  be  repeated  for  each  probe  position.  Once  hc  and  Tau>  have  been 
determined  and  either  p  or  pi  has  been  measured  at  the  same  probe  location,  then 
the  procedures  outlined  in  Section  3.1  can  be  used  to  obtain  the  flow  properties. 

Before  hc  or  Taw  can  be  determined,  the  thermal  accommodation  coefficient  of 
the  wire  must  be  known.  It  can  be  determined  from  Btatic  gas  tests  by  using  Eq. 
(19).  hc  is  determined  as  described  above  for  various  known  values  of  p  and  T  in  a 
static  vacuum  chamber.  For  a  100/i  inch  wire,  the  chamber  pressure  must  be  below 
1  torr  for  Eq.  (19)  to  apply.  Pressures  in  this  range  can  be  accurately  measured 
with  an  MKS  Baratron  pressure  transducer.  Measurements  at  several  values  of  p 
will  yield  an  accurate  measurement  of  the  effective  value  of  the  accommodation 
coefficient  of  the  probe  wire.  Unless  the  wire  is  used  to  measure  flows  that  will 
contaminate  the  wire  surface,  a  should  not  change  with  time.  However,  it  would 
be  helpful  to  occasionally  repeat  the  measurement  of  a. 

5.0  MEASUREMENT  RESULTS 


Preliminary  use  of  the  1/4  inch  free-molecule  hot  wire  probe  was  made  during 
the  plume  characterization  tests  in  the  4'  x  10'  vacuum  chamber  at  AEDC.  Mea¬ 
surements  were  made  in  the  plume  of  the  nozzle  (Fig.  2),  in  the  backflow  region, 
upstream  of  the  lip,  and  inside  the  nozzle.  The  measurements  were  made  without 
using  the  thermistors  as  heaters  to  eliminate  the  support  conduction  losses  nor  was 
the  probe  surrounded  by  a  radiation  shield.  Pitot  probe  measurements  were  made 
simultaneously  with  the  hot  wire  measurements.  Comparison  was  made  between 
these  measurements  and  measurements  of  the  number  density,  which  had  previ¬ 
ously  been  obtained  from  electron  beam  measurements  at  the  same  locations.  It 
was  demonstrated  that  the  hot  wire  measurements  were  able  to  resolve  the  structure 
of  the  flow  field,  although  the  magnitude  of  the  measurements  was  in  error.  The 
opportunity  did  not  exist  during  the  duration  of  the  tests  to  conclusively  demon¬ 
strate  the  ability  of  the  thermistor  heaters  and  radiation  shield  to  eliminate  the 
conduction  and  radiation  losses  and  additional  experiments  must  be  performed  to 
examine  their  performance. 

The  probe  was  operated  as  described  in  Section  4.3.  At  each  location  the 
probe  was  operated  at  5  or  6  overheat  ratios.  No  attempt  was  made  to  measure  the 
accommodation  coefficient  or  the  radiation  loss,  since  the  thermistor  heaters  were 
not  used.  The  reason  for  not  using  the  thermistors  or  a  radiation  shield  will  be 
discussed  at  the  end  of  this  section.  Because  a  and  Qj  were  unknown,  no  attempt 
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was  made  to  use  the  measured  values  of  hc  and  Tsw  to  determine  the  flow  properties 
but  rather  normalized  measurements  were  used  to  examine  the  flow  structure.  The 
measurements  were  also  valuable  because  they  provided  an  insight  into  the  operation 
of  the  probe  under  severe  temperature  conditions. 

The  flow  conditions  can  be  summarized  as  follows:  stagnation  temperature 
487°  to  500°K,  exit  Mach  number  4  and  tunnel  wall  temperature  (external  to  noz¬ 
zle)  77°K.  The  probe  operating  range  can  be  summarized  as  follows:  heating  power 
levels  4.9  x  10-#  to  1.432  x  10-4  watts,  average  wire  temperature  as  determined 
from  the  average  resistance,  324° K  to  562° K,  overheat  range  1.025  (for  measure¬ 
ments  inside  the  nozzle)  to  1.20  (for  measurements  in  the  backflow  region),  wire 
end  temperatures  measured  by  the  thermocouples  303°K  to  441°K,  probe  wire  adi¬ 
abatic  wall  temperature  within  the  isentropic  nozzle  core  flow  575.5°K,  measured 
adiabatic  wall  temperature  308.5°  to  543.7°K  and  measured  wall  heat  transfer  rate 
12.4  to  170.0  watts/m2oK.  The  probe  calibration  curve  is  shown  in  Fig.  20.  The 
probe  operating  conditions  and  measurements  of  Taw  and  hc  at  each  location  are 
given  in  Tables  I  and  II.  Typical  data  analysis  plots,  from  which  these  quantities 
were  determined,  are  shown  in  Figs.  32  and  33.  The  low  overheat  ratios  used  for 
some  tests  indicate  the  importance  of  accurately  knowing  i?2o  and  o2o  for  each  wire. 

The  measurements  indicate  that  large  temperature  gradients  existed  along  the 
wire.  These  gradients  resulted  in  large  conduction  and  radiation  losses.  The  extreme 
cases  occurred  for  measurements  that  were  performed  within  the  nozzle,  where 
the  difference  between  the  average  wire  temperature  and  the  end  thermocouple 
temperature  was  as  high  as  122“C  (run  number  40).  At  the  other  extreme  were 
the  measurements  performed  outside  the  plume  and  in  the-  backflow  region,  where 
this  temperature  difference  was  85°C  (run  number  28).  Both  of  these  temperature 
differences  occurred  at  the  maximum  overheat  ratio.  Clearly,  the  heat  losses  are 
much  too  great  to  expect  an  accurate  measurement  of  hc  and  Taw.  Therefore, 
thermistor  heating  will  be  required  to  eliminate  the  conduction  losses  during  future 
use  of  the  probe. 

Radiation  heat  transfer  from  the  surface  of  the  wire  and  the  black  resin  sur¬ 
rounding  the  thermistors  and  thermocouples  was  also  significant  under  the  test 
conditions.  This  subject  was  discussed  in  Section  4.0.  The  radiation  loss  from  the 
1/4  inch  wire  could  be  accurately  calculated  from  Eq.  (53)  if  the  wire  was  at  a 
uniform  temperature.  However,  since  that  was  not  the  case,  this  formula  could  only 
give  an  estimate  of  the  radiation  loss  by  using  the  mean  wire  temperature  as  deter¬ 
mined  from  the  mean  wire  resistance.  For  a  worst  case  estimate,  was  assumed 
to  be  77°K.  This  indicated  that  Qj  from  the  wire  was  in  the  range  of  3.2  x  10~6  to 
2.9  x  10-5  watts;  these  estimates  are  undoubtedly  high. 

Radiation  heat  transfer  from  the  surface  of  the  resin  blobs,  located  at  the  ends 
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of  the  wire,  was  more  difficult  to  estimate.  The  temperature  internal  to  the  resin 
was  known  from  the  thermocouple  temperature  but  the  surface  temperature  was 
unknown.  Assuming  that  the  resin  blobs  had  a  diameter  of  0.020  inch,  that  their 
surface  temperature  was  the  same  as  the  temperature  indicated  by  the  thermo¬ 
couples  and  that  the  surface  emissivity  was  1.0,  then  the  estimated  radiation  rate 
ranged  from  7  x  10-4  to  3  x  10“3  watts.  However,  this  is  considerably  greater 
than  the  power  supplied  to  the  probe,  indicating  that  the  surface  temperature  was 
much  less  than  the  indicated  thermocouple  temperature  and/or  the  emissivity  was 
much  less  than  1,0.  Also,  no  consideration  has  been  given  in  these  estimates  of  the 
radiation  shape  factor  between  the  walls  at  the  various  temperatures  which  existed 
within  the  4'  x  10'  chamber.  Nevertheless,  the  estimates  indicate  that  radiation,  as 
well  as  support  conduction,  produced  a  serious  heat  loss  from  the  probe. 

The  measured  profiles  of  hc  and  Taw  in  several  radial  planes  are  given  in  Figs.  34 
and  35.  Measurements  were  made  at  three  downstream  axial  stations  (5.74  inches, 
3.74  inches  and  0.74  inches,  respectively).  Some  measurements  were  made  inside  the 
nozzle  (-  0.26  inches)  and  in  the  backflow  region  upstream  of  the  lip  (-  0.26  inches). 
The  hot  wire  probe  was  mounted  adjacent  to  the  pitot  probe,  as  shown  in  Table 
II.  The  profile  at  5.74  inches  was  taken  across  the  entire  plume  profile  but  since 
it  was  symmetric  the  remaining  profiles  were  taken  only  from  the  jet  centerline 
(radial  position  3.48  inches)  outward.  The  measured  values  of  hc  and  Taw  were 
constant  through  the  core  of  the  plume  and  then  fell  rapidly  in  the  shear  layer.  The 
measurements  at  location  0.74  inches  indicate  that  the  isentropic  core  was  3  inches 
in  diameter  at  that  location,  which  means  that  the  shear  layer  was  0.5  inches  thick. 
This  is  approximately  the  thickness  of  the  boundary  layer  at  the  nozzle  exit.  The 
isentropic  core  then  diminished  and  the  shear  layer  greatly  thickened  as  the  flow 
expanded  downstream.  The  actual  flow  adiabatic  wall  temperature  was  estimated 
to  be  roughly  578°K.  This  value  was  only  approached  as  the  measurement  plane 
was  moved  upstream  and  inside  the  nozzle  the  measured  value  was  544°  K.  This 
lower  value  is  the  result  of  the  conduction  and  radiation  losses  from  the  probe.  In 
the  backflow  region  at  —0.26  inches  (runs  28  and  29),  hc  was  measured  to  be  12.9 
and  13.4  watts/m2oK,  respectively.  This  was  approximately  the  lower  limit  to  the 
measurable  heat  transfer  coefficient  and  resulted  primarily  from  radiation  rather 
than  convection  heat  transfer.  This  conclusion  is  also  confirmed  by  measurements 
made  from  a  dummy  probe  that  was  placed  in  a  vacuum  (p  =  4  x  10-3  torr) 
without  any  flow.  The  probe  had  not  been  calibrated  so  that  its  resistance, 
and  surface  area  were  not  accurately  known.  Measured  values  of  hc  with  this  probe 
in  the  vacuum  were  11.88  and  11.45/watts/m2oK.  These  values  are  very  close  to 
the  values  measured  in  runs  28  and  29. 

The  opportunity  did  not  exist  during  the  tests  to  fully  test  the  thermistor 
heaters  on  the  radiation  shield.  Thermistor  heating  tests  were  performed  on  an 
early  prototype  probe  which  had  a  1/4-inch  long  wire  of  diameter  2  x  10-4  inches. 
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This  probe  was  uncalibrated  and  Riq  was  only  approximately  known  (approximately 
50  ohm).  The  tests  demonstrated  that  the  thermistor  temperature  controller  could 
successfully  control  the  thermistor  power  to  place  the  thermocouple  temperature 
at  any  desired  set  point.  However,  when  power  was  applied  to  place  the  ends  of 
the  wire  at  the  average  wire  temperature  (as  determined  by  the  initial  average  wire 
resistance)  the  wire  power,  Qj,  was  reduced  to  essentially  zero  at  the  bridge  bal¬ 
ance.  In  other  words,  the  thermistors  were  heating  the  wire  as  well  as  the  wire 
supports.  Later  analysis  of  these  tests  indicated  that  this  was  caused  by  an  insuffi¬ 
ciently  accurate  knowledge  of  R20.  The  bridge  electronics  were  designed  for  a  probe 
having  a  resistance  of  10G0ofl  (1/4  inch  long,  50  x  10-6  inch  diameter  wire).  For 
such  a  wire  the  bridge  balance  resolution  would  be  0.05%.  For  the  200  x  10“  3  inch 
wire,  the  bridge  balance  resolution  was  only  1%.  Thus,  both  the  lack  of  knowledge 
of  and  the  resolution  of  the  bridge  balance  contributed  to  the  uncertainty  in 
the  average  wire  temperature.  In  the  tests  under  discussion  it  was  estimated  that 
the  thermistors  were  heated  to  7°C  above  the  average  wire  temperature.  This  is 
not  much  greater  than  the  expected  accuracy  of  the  bridge  electronics  (see  Section 
4.2.2).  Therefore,  the  thermistor  heaters  will  only  be  successful  if  the  temperature 
can  be  measured  with  great  accuracy,  and  if  the  resistance  bridge  has  great  res¬ 
olution.  This  will  require  using  a  small  diameter,  calibrated  probe  if  the  present 
electronics  are  used. 

The  radiation  shield  shown  in  Fig.  22  was  never  constructed.  However,  a 
preliminary  radiation  shield  test  was  performed  on  a  damaged  1/8-inch  probe.  The 
probe  was  surrounded  by  a  loosely-coiled  Nichrome  wire  helix  (0.03  inch  diameter 
wire).  The  helix  was  about  5/8  inch  in  diameter  and  2  1/4  inches  long.  The  shield 
proceeded  slightly  forward  of  the  wire.  The  probe  was  placed  downstream  of  the 
exit  of  a  low-density  sonic  orifice  {x/D  «  30,  Mach  number  at  probe  «  12  to  15).  A 
chr omel-  alumel  thermocouple  was  silver  soldered  to  the  Nichrome  wire  to  determine 
the  shield  temperature. 

The  relevant  parameters  of  the  test  are  as  follows.  The  gas  was  CO2  with  a 
stagnation  temperature  of  453°  to  463°K.  The  shield  could  be  heated  only  to  444°K 
with  the  DC  power  supply  available.  Without  thermistor  power  the  thermocouples 
on  the  wire  posts  measured  341°K.  With  maximum  thermistor  power  (12  x  10-3 
watts)  this  temperature  was  raised  to  390°K,  which  is  considerably  below  the  stag¬ 
nation  temperature.  However,  without  radiation  shield  heating  or  thermistor  power 
the  thermocouple  temperature  was  reduced  to  256°  K.  Therefore,  the  shield  raised 
the  thermocouple  temperature  85°  K  even  without  thermistor  power.  This  crude 
radiation  shield  had  a  noticeable  influence  on  the  temperature  of  the  ends  of  the 
probe  but  was  not  sufficient  to  eliminate  the  radiation  loss. 
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6.0  CONCLUDING  REMARKS 

The  development  of  a  free-molecule  hot  wire  probe  of  unique  design  has  been 
described  in  this  report.  When  coupled  with  an  independent  measurement  of  num¬ 
ber  density  or  pitot  pressure  it  is  capable  of  determining  all  of  the  flow  properties 
at  a  point  in  a  rarefied  flow  field.  It  has  good  spatial  resolution  and  large  dy¬ 
namic  range,  being  capable  of  obtaining  measurements  at  densities  below  which 
some  techniques  such  as  the  electron  beam  cannot  be  used.  It  can  also  be  used  to 
measure  flow  fields  within  nozzles  which  are  not  accessible  to  optical  and  electron 
beam  techniques. 

The  convective  heat  transfer  coefficient,  which  is  measured  by  the  probe,  is 
masked  by  influence  of  the  conduction  to  the  wire  support  posts.  Thermistor  heaters 
were  placed  at  the  ends  of  the  wire  to  add  sufficient  energy  to  just  compensate 
for  this  loss.  The  large  surface  area  of  the  resin,  which  held  the  thermistors  and 
accompanying  thermocouples  to  the  wire,  created  a  large  surface  which  could  lose 
energy  by  radiation  to  the  liquid  nitrogen  cooled  surfaces  which  surrounded  the 
flow  field  in  the  present  tests.  This  loss  must  also  be  eliminated.  That  can  be 
accomplished  by  using  smaller  thermistors  (0.005  inch  diameter  as  opposed  to  the 
0.013  inch  diameter  thermistors  used  for  the  present  tests)  to  reduce  the  radiating 
surface  area  and  by  using  more  efficient  radiation  shields.  Future  examination  of 
the  thermistors  and  radiation  shield  could  more  economically  be  performed  in  the 
free  -  molecule  flow  from  a  slit.  The  properties  of  this  flow  are  well-known  and  are 
given  in  the  Appendix  for  convenience. 

Two  electronic  control  circuits,  one  for  controlling  the  wire  temperature  and 
the  second  for  controlling  the  thermistor  power,  were  designed  and  constructed. 
The  tests  indicate  that  the  circuits  performed  very  well,  providing  the  accuracy  and 
control  required  for  the  proposed  use  of  the  probe. 

The  preliminary  tests  that  were  performed  with  the  probe  indicated  the  useful¬ 
ness  and  convenience  of  the  probe  for  making  flow  field  measurements.  Normalized 
measurements  compared  very  favorably  with  those  made  with  the  pitot  probe  and 
electron  beam.  Continued  development  will  be  required  before  the  probe  can  display 
its  maximum  utility. 
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P0  =  lOtorr 

Mq=4 

n^=  5.85  x  1021/m 

T0  =  500  K 

\=119K 

n  waii  =  1 0  /m 

Gas:  Carbon  Dioxide 

Pq_=  .063torr 

X<i  =  0.020  cm 

Figure  2. 

Nozzle  test  conditions. 

MODIFIED  STANTON  NUMBER 


Figure  3.  Modified  Stanton  number  for  a  cylinder  transverse  to 
the  flow  direction  in  free  molecular  flow. 
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Figure  9.  Relationship  used  by  Gottesdiener  (Ref.  2)  to  determine 
velocity. 
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Figure  13.  Detail  of  connections  at  end  of  wire. 


1/8  INCH  FREE  MOLECULE  HOT-VIRE  PROBE 
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Figure  18.  Measured  resistance  of  1/8  inch  probe. 
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Figure  21.  Measured  resistance  of  thermistors  on  1/4  inch  probe 
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Figure  24.  Hot  wire  control. 
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10k  Ik  100  10 

IL-0388 

NOTES: 

1.  The  decade  switch  Is  constructed  from  four,  10  digit  thumbwheel  switches. 

2.  The  numbers  to  the  left  ol  the  resistors  Indicate  the  valve  that  R  w  assumes  when  the  switch  is  in  that  position. 

3.  All  bridge  resistors  are  metal  film,  + 1%,  +  50  ppm/  ‘C,  1/10  wait. 


Figure  25.  Bridge  decade  resistor. 
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Figure  27.  Photograph  of  free- molecule  heat  transfer 
probe  electronics. 
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28.  Photograph  of  hot  wire  control  module. 
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Figure  29.  Photograph  of  proportional  temperature 
controller  module. 
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Measured  Wire  Resistance,  Rm  ohm 
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Input  Resistance,  R,  ohm 


Figure  30.  Bridge  calibration. 
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Input  Resistance,  R,  ohm 


Figure  31.  Bridge  calibration. 


Test  Run  26 

Overheat  Ratio,  1.19 

Heat  Transfer  Coefficient,  13. 3t o/m2K 

Adiabatic  Wall  Temperature,  324. 8 if 
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P ii  =  10.8  torr,  T„  =  487ff 
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Figure  32.  Probe  measurements  for  probe  outside  plume  shear  layer. 
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Heat  -  Transfer  Coefficient,  w/m2/K 
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Axial  Position 
Downstream  of 
Exit  Plane,  in. 


Radial  Location,  inches 

CG-O780 


Figure  34.  Radial  profiles  of  heat  transfer  coefficient  measured 
by  hot  wire  probe. 
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Axial  Position 
Downstream  of 
Exit  Plane,  in. 


Radial  Location,  inches 

CG-07S0 


Figure  35.  Radial  profiles  of  adiabatic  wall  temperature 
measured  by  hot  wire  probe. 
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Table  I 

Operating  Parameters  of  Hot  Wire 

Run  Number  Lowest  Overheat  Highest  Overheat 


Wire 

Heating 

Wire 

Heating 

Resistance 

Power 

Resistance 

Power 

(10_6)watts) 

(") 

(10_6(watts) 

Run# 

/Zi(Iow) 

AM 

A(high) 

A(M 

01 

271.0 

5.8 

286.1 

099.1 

02 

277.1 

5.8 

288.5 

096.8 

03 

281.4 

5.9 

290.9 

098.7 

04 

284.8 

6.0 

293.3 

107.4 

05 

285.7 

6.0 

293.3 

103.2 

06 

285.4 

6.0 

295.7 

126.3 

07 

286.2 

6.0 

294.7 

114.6 

08 

285.7 

6.0 

291.8 

086.2 

09 

283.3 

6.0 

290.9 

089.4 

10 

279.5 

5.9 

289.9 

097.0 

11 

274.8 

5.9 

288.5 

102.2 

12 

261.0 

5.6 

271.8 

055.1 

13 

251.0 

5.5 

259.8 

035.7 

14 

274.8 

5.9 

286.1 

081.5 

15 

282.9 

6.0 

293.3 

112.3 

16 

287.6 

6.0 

294.7 

L07.2 

17 

289.0 

6.1 

297.6 

137.2 

18 

288.0 

6.1 

297.6 

135.9 

19 

288.6 

6.0 

297.6 

135.8 

20 

288.6 

6.0 

297.6 

135.8 

21 

294.3 

6.2 

299.5 

139.4 

22 

294.3 

6.2 

299.5 

131.3 

23 

294.3 

6.1 

299.5 

131.1 

24 

294.8 

6.3 

299.0 

118.6 

25 

276.4 

5.9 

288.5 

088.5 

26 

223.8 

4.9 

245.5 

047.5 

27 

222.4 

4.9 

245.5 

047.6 

28 

22  L. 2 

4.9 

243.1 

047.7 

29 

221.2 

4.9 

247.9 

055.9 

30 

227.6 

5.0 

2503 

049.7 

31 

242.9 

5.3 

262.2 

057.3 

32 

272.1 

5.8 

288.5 

101.1 

33 

288.6 

6.1 

298.0 

117.9 

34 

295.2 

6.1 

300.4 

134.3 

35 

295.2 

6.2 

300.4 

134.9 

36 

295.2 

6.3 

300.4 

128.8 

37 

295.7 

6.2 

300.9 

131.3 

38 

297.6 

6.2 

302.3 

140.2 

39 

298.1 

6.3 

302.3 

138.6 

40 

298.1 

6.2 

302.8 

143.2  . 

41 

298.1 

6.3 

302.8 

141.6 
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Table  II 

Hot  Wire  Measurements  From  Nozzle  Tests 


Run 

Thermocouple 

hc 

Taw 

Location 

Number 

Temperature 

(watt/m2oK) 

(°C) 

Axial 

Radial 

(°c) 

(inch) 

(inch) 

01 

089.8 

046.3 

193.5 

5.735 

5.850 

02 

098.2 

054.0 

207.8 

5.380 

03 

106.4 

065.1 

221.2 

4.910 

04 

113.5 

078.2 

231.4 

4.440 

05 

117.0 

084.4 

234.6 

3.960 

06 

117.0 

088.0 

234.1 

3.490 

07 

118.1 

083.8 

235.6 

3.020 

08 

117.3 

083.7 

234.3 

2.540 

09 

111.6 

078.9 

229.9 

2.070 

10 

104.0 

057.6 

215.1 

1.600 

11 

097.4 

052.7 

201.9 

1.130 

12 

081.9 

031.4 

157.0 

0.183 

13 

073.8 

024.7 

132.4 

3.735 

0.183 

14 

101.6 

045.3 

201.3 

1.130 

15 

114.1 

069.5 

226.2 

1.600 

16 

125.4 

091.5 

239.9 

2.070 

17 

129.6 

098.4 

243.5 

2.540 

18 

128.7 

095.1 

242.8 

3.020 

19 

129.2 

095.1 

242.8 

3.490 

20 

130.1 

095.1 

242.8 

3.960 

21 

152.9 

155.8 

260.0 

0.735 

3.960 

22 

153.0 

141.4 

259.1 

3.490 

23 

153.4 

138.6 

258.7 

2.540 

24 

154.6 

141.7 

259.9 

2.070 

25 

125.0 

044.9 

206.3 

1.600 

26 

036.9 

014.7 

051.6 

0.145 

27 

035.5 

012.4 

043.1 

0.145 

28 

029.9 

012.9 

035.5 

-0.260 

0.145 

29 

030.0 

013.4 

039.6 

0.620 

30 

046.4 

014.4 

061.8 

0.735 

0.620 

'  31 

068.3 

019.8 

107.9 

1.090 

32 

105.7 

042.4 

196.7 

1.560 

33 

138.1 

080.1 

243.4 

1.800 

34 

151.5 

144.0 

261.5 

2.040 

35 

152.6 

144.5 

261.6 

2.270 

36 

152.9 

139.5 

261.9 

2.740 

37 

153.7 

136.6 

262.5 

3.220 

38 

165.9 

166.8 

269.4 

-0.260 

3.220 

39 

166.5 

170.0 

269.9 

2.900 

40 

166.8 

167.9 

270.5 

2.430 

41 

167.9 

164.5 

270.7 

1.900 
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APPENDIX  A 

FREE-MOLECULE  FLOW  PROPERTIES  DOWNSTREAM  of  a  RECTANGU¬ 
LAR  ORIFICE 

The  free-molecule  heat  transfer  probe  cannot  be  used  to  measure  a  flow  field 
without  an  independent  measurement  of  one  flow  variable,  such  as  the  density. 
However,  if  the  density  was  known  from  theory  for  some  flow,  then  the  operation  of 
the  probe  and  the  effectiveness  of  the  radiation  shields  could  be  verified  without  the 
need  for  a  companion  instrument.  The  free-molecule  flow  through  an  orifice  into  a 
vacuum  offers  such  a  flow  because  all  of  the  flow  properties  at  any  point  downstream 
of  the  orifice  are  known  from  free-molecule  flow  theory.  The  disadvantages  of  using 
such  a  flow  for  calibration  are  that  the  density  is  very  low  to  maintain  free-molecule 
flow  and  the  speed  ratio  is  much  lower  than  would  occur  in  the  case  of  a  nozzle. 
However,  it  does  offer  many  advantages  for  examining  the  probe  operation  which 
outweigh  these  disadvantages. 

Since  the  probe  is  a  long  cylinder,  a  rectangular  opening  was  chosen  for  the  ori¬ 
fice  geometry.  This  shape  will  result  in  flow  properties  which  are  reasonably  uniform 
over  the  length  of  the  wire  and  thus  would  require  less  correction  for  flow  gradients. 
While  the  properties  downstream  of  a  circular  orifice  are  well-documented,  those 
for  a  rectangular  orifice  could  not  be  found  and,  thus,  are  given  in  this  appendix 
for  later  use. 

We  will  consider  the  free-molecule  flow  from  a  gas  source  into  a  vacuum  (Fig. 
A-l).  The  orifice  has  dimensions  a  and  b.  Inside  the  source  the  velocity  distribution 

function  is  the  number  density  is  n„  and  the  temperature  is  Tv.  Let  £  be  the 
molecular  velocity  with  u  the  mean  velocity.  Then  the  number  density  at  any  point 
B  downstream  of  the  orifice  is  given  by  the  integral 

(A-l) 

where  the  velocity  distribution  is  given  by 

exphfe){Wvd*'£  ■  (A'2) 

The  Bolid  angle  subtended  by  dA  at  B  is 


duj  =  sin  tpd  <pd  9  = 


dAcosfi 

d* 


where 

d  =  a/  x'2  +  y'2  +  z'2  =  y/x2  +  y2  +  z2 


(A-3) 

(A-4) 
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cos  n  ~  —  . 

a 


(A- 5) 


Then  the  number  is  given  by 

"”(* e^{~Sr)^S^iA  (A‘6) 

A 

or,  losing  Eq.  (A-2), 


n  1  f  cosfi 

n„  4 tt  J  (P 


dA  . 


(A-7) 


The  geometric  integral  in  Eq.  (A-7)  was  evaluated  by  introducing  the  following 
coordinate  transformation  and  non-dimensionalization. 


x  =  x 

y\  =  ~(y  +  a/2)  y'2  =  ~(y  -  «/2) 

z[  =  -(z+b/2)  z’2  =  -(z-  6/2) 


X=*  Y=y-  Z=Z- 

a  a  a 

„  b  n  d 

R  =  -  D  =  - 

a  a 

D2  =  X2  +  Y2+Z2  . 


(A-8) 


(A-9) 


—  tan 


—  tan 


+  tan 


(A-10) 


The  result  of  the  integration  is 

i.±P(E+Jl^{*+r+*+^ ♦*>}-*) 

-  [n>-Y+nz  +  I(X  +  tf,p) 

(g,+y_flZ+l(1+^|-/zj 

-i  ^  ~  aX^  ~  R!2)  |  D2-Y-RZ  +  ^(1  +  JZ3)|  . 

The  x-component  of  the  velocity  at  any  point  B  can  be  calculated  from  the  integral 
nux  =  J f*  /„  df  =  £  cos  ft  (A-ll) 
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Using  Eq.  (A-2)  this  becomes 

“*  ("arav)  *  /  ^ <M  <A'12> 


(2irRTv 


\/2RT, 


Notice  that  the  velocity  relation  has  been  non-dimensionalized  in  terms  of  the  speed 
ratio 


5  = 


>/2M; 

The  evaluation  of  the  geometric  integral  gives  the  following  relation  for  the  velocity 
(or  speed  ratio)  at  any  point  B. 


4*^  ( —)  ~^=  = 
\rtvJ  y/2 BTV 


X2 

Z-R/2 

(Y-\)^X*  +  Y*- 

[  \^2+r2-r  +  i 

—  tan-1  | 

2  +  .R/2  \ 

^X*  +  Y2-Y  +  \) 

X 2 

tan  1  ^  — p== 

Z-Rf2 

yjX*  +  Y2  +  Y  +  i 


—  tan' 


Z  +  R/2 


-) 


,Jx2+Y2+Y 

(Z  -  aft)  r,„-i  / _ I+i _ ) 

VJr!+Z4--RZ  +  -R!/4)  [  \</X*+Z*-RZ  +  lP/iJ 


—  tan' 


'( 


y  _  I 

j* _ 2 


V^  +  Z2- J2Z-f-i?2/4 


)] 


(A- 14) 
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/X*  +  Z*  +  RZ  +  E?/  4 


w+m  L.-i  (  r+*  \ 

v'jc»  +  z=  +  flz  +  jp/4  [  \^/WTz^¥RzTWJl) 

f  .  r--  )1 

\y/x^rz*TW+Wji)\ 

+x (y+i  ~ y-$)  [s,n  (vJP  +  ^-flz+jp^) 


—  tan' 


sin  1  ^- 


/JT2  +  22  +J?2  +  fZ74 


K-  - 

*■  1  2  - 


+1/1  + 


^+i 


._  _i  (  I  x2  +  y*  -  y  +  i  V 
Sm  ^V2)2-F  +  i22  +  i(l  +  B)2J 

-■->  (J  JP  +  y  +  rtQl 

^Vj^  +  r+jjz  +  Jfi  +  jp)^ 

The  temperature  at  point  i?  can  be  determined  from  the  following  integral 

nRT  =  \j  (f -“)  /.<'?=  5  (A-15) 

The  first  term  can  be  evaluated  as  follows 

yV/v<*if=  (2^o)3/2  l  ?**P  —  dA 

A 

-  Sn'sr'  /  =  swnv 

/! 

when  use  is  made  of  Eq.  (A-7).  Then  the  temperature  distribution  can  be  written 


The  first  term  can  be  evaluated  as  follows 


nRT  =  nRTv  -  -nu^ 
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T 

T„ 

T_ 

Tv 


=  1 

=  1 


i  ji. 

3  RTV 


(A-16) 

(A-17) 


The  number  density  profiles  in  the  plane  Y  =  0  are  shown  in  Fig.  A-2  for  four 
distances  downstream  of  the  orifice  plane.  The  number  density  has  been  normalized 
by  its  value  upstream  of  the  orifice.  The  profiles  are  reasonably  flat  over  the  length 
of  the  wire  (to  Z  —  5.2).  Fig.  A-3  shows  the  corresponding  speed  ratio  profiles  at 
the  same  locations.  Notice  that  the  speed  ratio  is  always  less  than  1.0;  at  large  X 
it  would  approach  a  limit  of  2/v'jF  on  the  centerline.  Therefore,  the  speed  ratio  can 
never  approach  the  high  values  observed  in  the  nozzle  find  plume  and  this  is  one  of 
the  major  disadvantages  of  using  this  flow  for  calibration  puiposes.  The  advantages 
are  that  the  density  is  well  known  at  any  location  from  the  theory  so  that  the  probe 
operation  can  be  verified  without  a  separate  density  measurement. 
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z 


Figure  A-2.  Number  Density  Profiles  Downstream  of  a  Rectangular 
Orifice  for  Free  •  Molecule  Flow. 
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Figure  A— 3.  Speed  Ratio  Profiles  Downstream  of  a  Rectangular  Orifice 

for  Free  -  Molecule  Flow. 
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NOMENCLATURE 

Function  defined  in  Eq.  (33) 

Exit  area  of  nozzle 
Surface  Area  of  wire,  irDL 
Throat  area  of  nozzle 
Function  defined  in  Eq.  (34) 

FUnction  defined  in  Eq.  (39),  A  +  <xrB 

Isentropic  thrust  coefficient 

Measured  thrust  coefficient 

Specific  heat  at  constant  volume  of  gas  at 

static  temperature 

Specific  heat  at  constant  pressure  of  gas  at 
static  temperature 

Average  molecular  speed,  yj ^RT 

Wire  diameter 

FUnction  defined  in  Eq.  (13) 

Function  defined  in  Eq.  (5) 

Function  defined  in  Eq.  (6) 

Enthalpy  of  gas  at  static  temperature 
Convective  heat  transfer  coefficient  from  wire, 

Qcf\As{Pw  ~~  Taui)] 

Wire  current 

Modified  Bessel  function  of  first  kind  and  order  v 
Number  of  kT/2  units  of  internal  energy  carried  by  molecule 
Thermal  conductivity  of  gas  at  static  temperature 
Function  defined  in  Eq.  (10) 

Knudsen  number,  X/D 

Thermal  conductivity  of  wire  material 

Boltzmann  constant 

Length  of  wire 

Fiee  stream  Mach  number 

Mass  of  gas  molecule 

Function  defined  in  Eq.  (44) 

Gas  number  density 
Nusselt  number,  kcDjK 
Static  gas  pressure 
Rocket  motor  chamber  pressure 
Prandtl  number  of  gas  (icp/K 
Rate  of  heat  transfer 
Gas  constant  of  gas 
Wire  resistance 
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Re 

r 

St 

s 

T 

T, 

V 
Vm 

x 

X 

Y 

Zi,Z2t 

a 

atr 

e 

e<v 

7 

A 

P 

P 

a 

9 


Reynolds  number,  pvD/fi 
T  —  T 

Recovery  factor,  — — 

T0 -T 

Stanton  number,  NuRejPr  =  hcfpcpV 
Speed  Ratio,  V((2 RT)1!2 
Temperature 

Temperature  at  junction  of  wire  and  support  posts 
Velocity 

Mean  molecular  speed,  \/2RT 
Coordinate  along  wire 
Rinction  defined  in  Eq.  (41) 

Function  defined  in  Eq.  (41) 

Functions  defined  in  Eq.  (3)  and  (4) 

Thermal  accommodation  coefficient  , 

Temperature  coefficient  of  resistivity  of  wire  material 
A  t/At 

Emissivity  of  wire  material 
Ratio  of  specific  heats,  cpfcv 
Mean  free  path 

Gas  viscosity  at  static  temperature 
Gas  density,  rnn 
Stefan-Boltzmann  constant 
Tw  -  Taw 


SUBSCRIPTS 


a 

aw 

c 

i 

i 

k 

o 

r 

t 

tr 

w 

«"oo 

OO 


Reference 

Adiabatic  wall  (equilibrium) 

Convection 

Impact 

Joulean  heating 

Conduction 

Stagnation 

Rotational 

Thermistor  heating 

Translational 

Wire 

Infinitely  long  wire 

Solid  Walls  surrounding  probe 
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